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. Abstract. A method for measuring the transmittivity of optical samples by using 

squeezed-vacuum radiation is illustrated. A squeezed vacuum field generated by a 
below-threshold optical parametric oscillator is propagated through a nondispersive 
, medium and detected by a homodyne apparatus. The variance of the detected 

' quadrature is used for measuring the transmittivity. With this method it is drastically 

reduced the number of photons passing through the sample during the measurement 
interval. The results of some tests are reported. 

+-> 

c 

O" 1 ' 1. Introduction 

Vacuum fluctuations of electro-magnetic (e.m.) fields have been the ultimate limit on 
the precision of optical measurements until the advent of squeezed light. Since then, 
many attempts have been made for reducing the shot noise level blurring several types 
of signals. Caves £Q first proposed to combine coherent and squeezed vacuum radiation 
for overcoming the quantum limit in gravitational wave antennas. Grangier et al. j2] 
up-graded a polarization interferometer by injecting a squeezed vacuum through a dark 
port. Polzik et al. [3] provided stunning evidence that a gain of some dB over the 
standard quantum limit is achieved in the resonant interaction of atoms with squeezed 
light. Their experiment was performed by combining in a well defined phase relation 
a coherent field with the output of an optical parametric oscillator (OPO) operating 
below-threshold. 

Other experiments have exploited the correlation between twin beams for reducing 
the noise level of the probe field [3J U\ ■ The principles of these measurements were 
highlighted in Ref. |3J, where it was recognized that the losses occurring in one beam 
can be inferred from those relative to the other one (see also jEllZllH])- 

In this communication a method for probing the transmittivity T of a sample with 
squeezed vacuum radiation is discussed. Standard methods rely on direct measurement 
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of the radiation intensity entering and leaving the sample. Sufficient accuracy can 
be achieved by using beams so intense to contrast the shot-noise, although, in some 
circumstances, using high input intensity is either not useful (in case of very low 
absorption) or unwise (strongly non-linear materials or samples whose structure may 
be altered by intense photon fluxes). 

In alternative to the above schemes PjJ E] in the proposed method the sample 
is irradiated with a squeezed vacuum field. Then, the emerging one is combined with 
a coherent one (local oscillator, LO) in a balanced homodyne detector measuring the 
fluctuations in a suitable spectral range. The interaction of the squeezed vacuum with 
the sample modifies the spectrum of the homodyne current by changing its variance. 
Hence, the transmittivity is determined by measuring the variance changes. The main 
advantage of this method is a very low number of photons interacting with the sample. 

Below-threshold degenerate OPOs produce e.m. radiation represented by a 
combination of squeezed vacuum and thermal components (squeezed-thermal-vacuum 
states, STV) with a Gaussian statistics. The OPO working conditions determine the 
STV state properties EI] • The propagation through non resonant media transforms 
such a state into another STV one with different variances of the field quadratures AX^. 
The change of AX| is used for measuring the transmittivity T (Section 2). 

The squeezed radiation is analyzed by a balanced homodyne detector providing 
the field quadratures X^ = | {ae~ l< ^ + a^e*^) via the controlled interference between the 
STV state and a strong coherent LO of relative phase (p. Since the detected signal is 
proportional to X^ times the LO amplitude, the detection is efficient also in case of 
very weak beams, as in the present case. Consequently, the effects of the SNR on the 
accuracy can be disregarded. 

Essential to this method is the use of a Gaussian distributed quadrature X^. 
This means that for testing the method it is necessary to preliminarily measure the 
distribution function by sampling X^ an adequate number N of times. In alternative, 
it is also possible to determine the whole Wigner function with quantum homodyne 
tomography (QHT) jTJ] using samples uniformly distributed over the whole interval 
(0, 2tc). Distributing N samples in the interval (0, 27r) reduces the accuracy of only a few 
percent. This slight loss is largely compensated by a three dimensional characterization 
of the STV state in the phase space. 

Aim of this communication is to assess the feasibility of this scheme by testing the 
validity of two main assumptions, namely i) the generation of Gaussian STV states by 
a below-threshold OPO, and ii) the description of the absorption process as a simple 
scaling of the P-representation of the STV state. The dependence of the accuracy of 
the proposed method on the STV state parameters is also examined. Moreover, the 
method accuracy is compared with that achievable with standard techniques. Some 
measurements carried out with a below-threshold type-I Lithium Niobate (LNB) OPO 
at A = 1064 nm, typically generating few pW STV states are illustrated. 

The paper is organized as follows. In the next section the properties of the 
STV states undergoing lossy propagation are discussed. Then, in Section 3, the 
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accuracy of the measurement of T based on this method is compared with that of 
different techniques. Section 4 is dedicated to the description of the experimental tests. 
Eventually, in Section 5, conclusions are drawn. 

2. Generation and propagation of STV states 
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State Quantum 
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Figure 1. Schematic of the OPO cavity and the experimental setup. The STV states, 
generated by the OPO, pass through a sample of variable transmittivity T and then 
are characterized by a homodyne detector. State parameters are recovered by QHT 
data processing. 



The field generated by a below-threshold OPO satisfies the Langevin equation: 



(k + iip + iSip (t)) a 
(k — iif> — i5if> (t)) at 

(1) 

where a™ is the noise entering the cavity through the output mirror M\ (see Fig. 1) 
and a % 2 represents both the noise entering through the input mirror M2 and the crystal 
contribution. The quantities K\, k 2 are damping coefficients, whereas k = K\ + k 2 . The 
parametric gain is the sum of a stationary mean value, 7, and a small fluctuating 
contribution, £7 (t) . Similarly, if} and Sif> (t) are the mean cavity detuning and its 
fluctuations, respectively. In this context, the ratio £ = 7 2 //t 2 represents the distance 
of the actual operating condition from the OPO threshold power, while K\/k is the so 
called coupling efficiency. 

In the frequency domain the variance AX 2 of the output quadrature X = X^ =0 
reads as: 
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where u is the frequency offset from the optical frequency Uq, and £7 (t) and 5ip (t) have 
been neglected. The variance AY 2 (Y = X^, =7r / 2 ) is given by a similar expression with 
7 replaced by —7. 

For a single-input cavity (k 2 = 0), and for ip = 0, the product 16AX 2 AF 2 reduces 
to unity, corresponding to a minimum uncertainty state. In general, this condition is no 
more satisfied for double-ended cavities or non-zero detuning or lossy crystals. Then, 
it is worth characterizing the OPO output at the sampled frequency u by means of the 
adimensional parameters: 

nth = 2 (VaX 2 AF 2 - ij 

1 / I AX 2 I AY 2 \ 
"* = 7 Va^ + Va^5- 2 > ( 3 ) 



representing the average number of thermal and squeezed photons, respectively. They 
measure the deviation of the actual state from the minimum uncertainty one and its 
effective squeezing. In particular, the mean total photon number is given by: 

N t ot = n sq + nth + 2n sq n th , (4) 
while the variance of the generic quadrature reads: 



AXl= ( 2nth + V ( 1 2;^ 2 V i] //,.„ ccs 2o ) . (7) 




0.6 0.7 0.8 0.9 0.6 0.7 0.8 0.9 



Kj/K Kj/k 

Figure 2. n sq (a) and nth (b) vs. coupling efficiency k\/k at half of the threshold 
power {£ = 0.5) for to = 0. The curves refer to three different detunings ip — 
0.0, 0.1, 0.2 (full, dashed and dot-dashed lines). 

The photon numbers n t h and n sq depend on the frequency offset u and on the 
OPO parameters, namely, distance from the threshold (S), escape efficiency (kx/k), 
and cavity detuning (ip). In Figs. 2-a and 2-b n sq and n t h as functions of the ratio 
K\/k and for three different detunings (ip = 0.0, 0.1, 0.2) have been plotted (OPO gain 
one-half of the threshold £ = 0.5 and u — 0). The detuning plays a more significant 
role in proximity of ki/k~ 1, that is for a single-ended cavity configuration. 

Since these states have been obtained by neglecting the time-dependent part of both 
gain and detuning, they share the Gaussian statistics of the driving fields a l ™ 2 . They can 
be imagined as obtained by squeezing a thermal state v = (nth + I) -1 [nth/ (nth + l)] at< \ 
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with n t h the number of thermal photons (see Eq.([3]-a)), whose actual temperature is 
not necessarily coincident with the local one. As a consequence the density matrix is: 



Q 



S{C)vS\Q , (6) 



where S(() = exp{^C a ^ 2 — |C* fl2 } is the squeezing operator [n sq = sinh 2 \(\ ) P ITU]. 

These STV states are described by a Gaussian Wigner function centered at the 
origin: 



W («) = " — T^r^T* ex P 



2vrVAX 2 Ay 2 "V 2AX 2 2AF 2 
- J P(/3)exp (—2 |« — /3| 2 ) d 2 p, 



with P (/?) the corresponding P-representation: 



P (/?) = _ exp 



3 



2vr v /(AX 2 -i) (AF 2 — |) V 2 (AX 2 - I) 2 (AF 2 - J) 
After propagation through a medium of transmittivity T the density matrix 
P (a) | a) (a| d 2 a 



modifies as 



g T = / P (a 



VTc^I (^VTa d 2 a = J Pt (a) |a) (a| d 2 a . 



with 



1 (a 
P T (a) = -P 



T \Jt 



3? [a] 2 3 [a] 



2^{AX T -\) (AF T 2 -I) 6XP V 2 (AX 2 -I) 2 (AF r 2 - J) 



and 



AX| - ~ = T (AX 2 - . 

AX 2 — | is the deviation of the actual STV variance from the vacuum state case (shot- 
noise). A similar expression is found for AY^. 

In principle, in the absence of multiple reflections within the sample, the 
transmittivity T is given by T = T{T s i a {T2^ where T\ and T2 are the Fresnel transmission 
coefficients at the input and output faces of the sample respectively and T s[ab is the 
sample internal transmittivity. 

Next, introducing the subfixes and T for labelling up- and down-stream 
quantities, respectively, for a generic quadrature X^ the variance transforms as: 



1 „/.„o 1 

4 



AXl T -- = T AXj - 7 ) . (7) 
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A 


B 




B(Q HT ) 


N tot 





1 


-0,05 ±0,07 


1,1 ±0,1 


n t h 


0.12 


0.89 


0.07 ±0.05 


0.85 ±0.07 


n S q 


-0.12 


1.14 


-0.16 ±0.05 


1.14 ±0.07 



Table 1. Coefficients A and B computed by Eas.l jl()ll2(l (left) and experimental ones 
measured by QHT (right). 



Accordingly, T can be obtained by measuring the up- and down-stream quadrature 
variances: 

This relation suggests a simple way to measure T through the deviations of a generic 
quadrature from the vacuum noise level. 

By means of Eqs.© and (jHJ), T can be also expressed as: 

(2n t h,T + 1) (1 + 2n sq , T + 2J(1 + n sq , T ) n S g iT cos20) - 1 

T= ) { • (9) 

(2n thfi + 1) ( 1 + 2n sqfi + 2^/(1 + n S g )0 ) ™ S g,o cos 20) - 1 

On the other hand, N tot transforms proportionally to T as for a classical field: 

N tot , T = TN totfi . (10) 

Using Eq.(JD) in the above expression and combining it with Eq.® n thjT and n sq ^ T can 
be expressed in terms of T and of the initial values n t h,o and n S(Ji o : 



2n thjT + 1 = W [1 - T ± T(l ± 2n <M ) (1 + 2n s ,, )] 2 



n 2 



2T(1 + 2n t h,o)\ (1 + rc sg ,o) ™ sg ,o 



1-T + T ( 2^ + !)(! + 2^ ) 
9 ' 2n^ T ±l 1 ^ 

For the STV state used in the test discussed in Section 4 (n t h,o = 0.55 and n sq $ = 0.11) 
nth,T and n sq ^ are practically linear in T, that is: 

^ = A,±^T 
nthfl 

^L = A sq + B sq T. (12) 

In Table ^ the coefficients calculated by linearizing Eqs. ()11|) have been reported (first 
two columns) together with those obtained experimentally (last two columns). For the 
sake of completeness the measured ratio N toti T / N tat>0 has been reported as well, in order 
to evidenziate the agreement with the theoretical value of Eq.([10|). 

In the measurements discussed below, T was determined through a direct 
measurement of the parameters n t h and n sq by QHT technique based on pattern 
functions. 
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In conclusion, it is worth remarking that the above expressions of T are valid 
for Gaussian field quadratures. A Gaussian statistics follows from the assumption of 
time independent gain and detuning of the OPO. In Section 4 the correctness of this 
assumption will be discussed for the used OPO by measuring the deviations from the 
Gaussian statistics by means of the kurtosis parameter if^, vanishing for the Gaussian 
case, defined as: 

AX* 



AXf being the fourth order moment of Xa 



3. Accuracy 



The limit of the uncertainty on the estimate of T expressed by Eq. ( 
confidence interval 5 [AX|] in the measurement of AX|: 



depends on the 



6T 



Since 



AX 



5 [AAJJ 2 r _ 2 



(14) 



with N the number of acquired data, the relative error on T is given by: 



5T 



N lAXjo 




+ 



AX lo 



1 

T 



(15) 

This expression gives, for a given T, the relative error as a function of N and AX^ which 
in turns depends on the OPO working condition, namely, distance from the threshold 
(£), escape efficiency (ki/k), and cavity detuning (-0). On the other hand, the total 
number of photons N p h hitting the sample during the measurement is: 



N 



ph 



N tot NKT x 



(16) 



with t s 1 the sampling rate and N tot given by Eq.(jH). 

In Fig. 3 Eq. I15p ) has been plotted as a function of the transmittivity 

T for uj = ip = 0, £ = 0.5, and three different escape efficiencies (k%/ k = 0.5, 0.75, 
1). The relative error increases for T approaching zero. Fig.3-b gives the photon dose 
Nph (Eq. (|16|)) necessary to obtain a relative error ^ = 0.01 for the parameters of 
Fig.3-a and kt s = 6. The plot evidentiates the increase of N p h by more than an order 
of magnitude for T less than 0.01. 

Instead of keeping constant during the acquisition of the N samples, the angle 
can be varied uniformly in the interval < <fi < 2tt. The ensemble so obtained 
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Figure 3. (a) Relative error f fEci.l|15|l) and (b) number N p h of photons hitting 

the sample for ^ = 0.01 (Eq.^SJ and kt s = 6) vs. transmittivity T. The plots refer 
to u) = ip = 0, at half the threshold {£ = 0.5) and ki/k = 0.5, 0.75, 1) (dot-dashed, 
dashed and full lines). 



can be processed by means of QHT ^T] for obtaining the field Wigner function. 
The tomographic processing can be based on the so-called pattern function method, 
consisting in averaging the pattern function R v [0] (Xg, 9j) relative to an assigned 
operator O, and having for argument the j-th realization Xg. of Xg for the LO phase 



9j, 



i N 



N , 



The subfix "rj" indicates the dependence of the pattern function on the homodyne 
efficiency rj. 

For the operator O = AX? the confidence interval provided by this method reads 



5qht [AX}} = -^V Ai?2 i AX l1 ' (" 



with AR 2 [0} = R 2 V [0] - R v [0] . Consequently Eq. dH is still valid with 5 [AX}] 
replaced by 5 QHT [AX}] . 

Next, taking into account the analytic expressions of i?^[AX|] and R V [AX}} [T2j, 
it can be shown that: 



AR 2 



AXi 



C + d cos (20) + C 2 cos (40) , (19) 



with the coefficients Co, C\ and Ci given in Appendix. For the variances AX 2 , AY 2 
relative to OPO devices similar to that used in the experimental test, 5qht [AX 2 ] 
differs from 5 [AX?] only by some percents. This means that collecting N samples in 
the interval (0, 2tt) reduces the accuracy with respect to the constant phase case by 
only a few percent. 

Conventional measurements of T using coherent CW probe beams and the radiation 
power, P, as observable, are in some way corrupted by the detector noise equivalent 
power (NEP), and the measurement error reads: 



5P = ^/fvuooB P + NEP , (20) 
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with u the radiation frequency, and B the detection bandwidth. 

Measuring T as the ratio Pt/Po of the power down- and up-stream the sample the 
relative error is: 



ST 
~T 



SNR 



\ 



J^2 



1 + 



hm B SNRT 
NEP W 



+ 1 + 



fkvoB SNR \ 
NEP N J ' 



(21) 



with SNR = P /NEP and N the number of data. 

The total number of photons (see Eq. ([16))) passing through the sample during the 
measurement interval is now given by 



N 



ph 



SN R ^Nr s 



(22) 



so that, the factor j^Sr* in Eq. (jUJ can be replaced by sn * Bt ° (with Bt s >1). Then 



ph 



the ratio is a function of T and SNR. Using for SNR the limiting value 



ph 



the plot of Fig. 4 representing N p h vs. T for 4jr = 0.01 and 5r s = 10 has been obtained. 
Comparing it with Fig. 3-b it appears evident the much lower photon dose required by 
the present method. 
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Figure 4. Photon dose passing through the sample {N p h) vs. T in case of standard 



transmission measurements. The curve refer to Bt s 



10 and = 0.01. 



In case a coherent beam (a a ) is mixed with a squeezed thermal vacuum one (asTv)-, 
as in Ref . jS] , the total field is described by: 

d-tot = e Qstv + a a , 

with 9 their locked phase difference. 

It can be shown that the addition of the squeezed component modifies slightly Eq. 
(120)) with B replaced by 

B eff = B ( 1 + n sq + n th + 2n sq n th + J (1 + n sg ) n sq cos 26 1 
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For COS26 1 = —1 and a/(1 + n sq ) n sq > n sq + n th + 2n sg n t h the squeezed vacuum 
component reduces the effective detector bandwidth. The reduction of B (typically 
B e ff > .5B) implies a proportional decrease of N p h for assigned %f and T. 

4. The experiment 

The reliability and accuracy of the method were tested with a sample of variable 
transmittivity. The T values obtained via QHT were compared to those measured, 
with an accuracy of 10~ 4 , with standard techniques employing 1 mW coherent beam at 
A = 1064 nm. A schematic of the experimental set-up is shown in Fig. 1. 

STV states were generated by a degenerate type-I OPO and characterized by a 
homodyne detector, both described in details in Ref. [T3|. In the present case, cavity 
mirrors were adjusted in such a way as to have a cavity linewidth of 15 MHz. 

The OPO output was propagated through a variable neutral density filter, which 
changes T without introducing misalignment, and keeping homodyne visibility at a 
constant value. The transmittivity T was varied between 0.45 and 1 in discrete steps. 
The beam passing through the non-absorbing zone (T — 1) of the filter was used as 
a reference state. The field leaving the absorber was sent to an homodyne detector 
with an overall efficiency of 77 = 0.88 ± 0.02. The average electrical signal level at the 
homodyne output was 15 dB higher than the electronic noise. 

Tomographic data were acquired by sampling the homodyne signal. To avoid any 
effect of the laser technical noise on the measurement, data sampling was performed by 
mixing the homodyne current with a sinusoidal signal of frequency Q = 5 MHz. Then, 
the resulting current was low-pass filtered, with a cut-off frequency of 2.5 MHz, and 
10 6 samples were collected with at 2.5 Msample/s (t s = 400 ns) in order to pick-up 
statistically independent data. 




-0.06 -0.04 -0.02 0.00 0.02 0.04 0.06 -0.10 -0.05 0.00 0.05 0.10 



Figure 5. Distribution of X values measured for £ = 0.5 a) and 0.95 b). The kurtosis 
Kq (Eq. dU) is respectively equal to 0.005 and 0.5. Full lines represent Gaussian with 
the same mean and variance. 

Fixing S = 0.50 the reference STV state had N totfl = 0.79±0.06, n thfi = 0.55±0.02 
and n sq fi = 0.11 ± 0.01, corresponding to a photon flux of 10 7 s -1 . For this state it 
resulted < 0.01 (see Fig. 5-a and Eq. (fTT?|) ) for any (ft, thus indicating that the 
corresponding quadrature statistics was very close to the Gaussian one. 
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In order to reduce the influence of residual fluctuations of the STV state, each 
experimental point was averaged over multiple (~5) tomographic acquisitions. In the 
present conditions the QHT error was negligible with respect to the standard deviations 
of the STV state parameters. 

To assess the robustness of the method, the transmittivity, Tqht, obtained by 
tomographic reconstruction was compared with the corresponding value, T st , provided 
by standard intensity measurements. 




0.4 



0.5 0.6 0.7 0.8 0.9 1.0 

T 

Figure 6. Tqht vs. T st . Experimental points are plotted together with the expected 
behaviour of Eq. (|10f> (straight line) . 



In Fig. 6, Tqht = N to t,T / 'N to t,o (see Eq. (fTTI|l ) was plotted vs. T st together with 
the expected behavior Tqht — T st (straight line). A linear regression of the data with 
Tqht = Af ot HT) + B$ HT) T st7 gave 4f T) = -0.05 ± 0.07 and B% HT) = 1.1 ± 0.1 in 
good agreement with the expected values of A tot = and B tot = 1 respectively. 

In order to estimate Tqht through other quantities, the measured value of 
n S q,T/n sq ^ versus T st was plotted in Fig. 7 together with the linear approximation 
of Eq. ()12l -b). Linear regression on experimental data gave A^ HT) = -0.16 ±0.05 and 
g(QHT) _ q.07, va l U es in good agreement with A sq = —0.12, B sq = 1.14. 

Each experimental point of Fig. 7 represents an average value obtained over 
multiple acquisitions. In the inset the different values of Tqht-, corresponding to four 
acquisitions at T st = 0.64 are reported. The bar indicates the quantum limit error, 
calculated by using Eq. (fT3j) . As it can be seen, all the points are spread over a range 
comparable to the quantum limit. 

In Fig. 8 the shadowed area represents the quantum limit for the accuracy vs. T 
(see Eq. (jl5|) ) for the present experimental conditions and N = 10 4 . The accuracy width 
is almost constant in the tested range of T while it deteriorates for low transmittivity, 
as expected. 

Finally, an identical behavior was observed for nth,T / nth,a ( n °t plotted) resulting in 
A {QHT) = o.o7 ± 0.05, B\2 HT) = 0.85 ± 0.07 (A th = 0.12, B th = 0.89). 
A summary of the experimental findings is reported in Table 
The photon flux at the OPO output F = N tot /r, with r the cavity photon lifetime, 
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Figure 7. n sq _T vs. T st . Experimental points are compared with Eq.H12b). The 
points in the inset are relative to four measurements for an attenuator transmittivity 
equal to 0.64 while the error bar has been calculated by using Eq. (fTSft . 
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Figure 8. Quantum limit (shadowed area) of the relative error on T for the 
experimental case discussed in the text and N = 10 4 . 



was less than 10 7 s -1 , for N to t ^ 0.7 and r pa 6.6 x 10 -8 , corresponding to an optical 
power < 4.2 pW. The method was tested for different input states, by varying S and 
hence the photon flux by showing a good reliability down to a photon flux F~5x 10 6 
s- 1 (i.e. ~2.2 pW and N totfi = 0.37). 

With N = 10 6 8 QHT [AX 2 ] ~ 1. 3x 1(T 3 and 5 QHT [AY 2 ] ~ 0.8 x 10^ 3 corresponding 
to 5T/T ~ 0.0024 and ~ 0.056 for T — 1. These QHT estimates were slightly less 
accurate than those one could obtain by concentrating N/2 data on X and N/2 on Y 
quadratures and computing their variances. 

5. Conclusions 

A scheme for measuring the optical transmittivity of a sample by using squeezed vacuum 
radiation has been illustrated. Main advantage of this method is a number of photons 
hitting the sample during the measurement some orders of magnitude smaller than that 
relative to standard techniques based on intensity measurements of coherent beams. 
The core of the method consists in the measurement of the variance AX? of a 
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generic quadrature of a squeezed vacuum field, generated by a below threshold 0P0 
and passing through the sample under investigation. The quadrature is measured by a 
homodyne detector. In the simplest implementation AX? is obtained by averaging the 
squared samples relative to a constant phase <fi. In the test described in the paper 
has been obtained by scanning the interval (p G (0, 2n). This approach has been 
preferred since it provides a complete reconstruction of the squeezed vacuum Wigner 
function. 

Essential to this scheme is the assumption of Gaussian statistics for the squeezed 
vacuum field. This property has been checked on the recorded samples relative to a 
given phase and confirmed by the field Wigner function. 

The accuracy of this method has been compared with that based on absorption 
of coherent beams (with and without a squeezed vacuum component) as a function 
of sample transmittivity, number of data and detection bandwidth. In the case the 
number of photons interacting with the sample during the measurement is an important 
parameter, the proposed method is the most accurate. 

The experimental tests have shown that, for photon fluxes of the order of few pW 
(at 1064 nm), the accuracy is of the order of the quantum limit, that is the method does 
not suffer substantially from other technical noise sources 
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Appendix 

The coefficients C , C\ and C 2 of Eq. (fT^j) are given by: 



C = y (AX 4 + AF 4 ) + 9AX 2 AF 2 + (\ 

d = - (AX 2 - AY 2 ) [3 (AX 2 + AY 2 ) - l] 
C 2 = l(AX 2 -AY 2 ) 2 . 



V 



3 



) 



(AX 2 + Ay 2 ) + 




) 



